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Abstract 

We present a determination of the isovector, P-wave nn scattering phase shift obtained by extrapolating recent lattice QCD re¬ 
sults from the Hadron Spectrum Collaboration using m n = 236 MeV. The finite volume spectra are described using extensions of 
Liischer’s method to determine the infinite volume Unitarized Chiral Perturbation Theory scattering amplitude. We exploit the pion 
mass dependence of this effective theory to obtain the scattering amplitude at m n = 140 MeV. The scattering phase shift is found 
to agree with experiment up to center of mass energies of 1.2 GeV. The analytic continuation of the scattering amplitude to the 
complex plane yields a p-resonance pole at E p = ^755(2)(1)(q°) - ^ 129(3)(l)(j)] MeV. The techniques presented illustrate a pos¬ 
sible pathway towards connecting lattice QCD observables of few-body, strongly interacting systems to experimentally accessible 
quantities. 
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The spectrum of hadronic resonances has long served as a 
window into the non-perturbative nature of Quantum Chromo¬ 
dynamics (QCD), the fundamental theory of the strong force. 
Hadronic resonances are color-singlet combinations of the fun¬ 
damental degrees of freedom of QCD (quarks, anti-quarks, and 
gluons). They are observed as unstable resonant enhancements 
in the scattering of QCD stable hadrons, such as the pion. A 
simple example of a hadronic resonance is the p that occurs in 
nn: scattering. The non-perturbative nature of QCD makes di¬ 
rect determination of the properties of hadronic resonances a 
challenging task. 

Presently, the only means to study properties of low-energy 
hadronic states in a systematically improvable way is to per¬ 
form a non-perturbative numerical evaluation of the QCD path- 
integral, by statistically sampling the gauge fields in a dis¬ 
cretized finite volume to obtain correlation functions. This pro¬ 
gram is known as lattice QCD. The last decade has witnessed a 
tremendous advance in the ability of the lattice QCD commu¬ 
nity to connect experimental phenomena directly to the stan¬ 
dard model of particle physics. It is not unreasonable to expect 
that in the upcoming decade most “ simple ” observables, such as 
masses, decay constants and elastic form factors of low-lying 
QCD stable particles, will be computed using physical values 
of the quark masses and QCD+QED gauge configurations (see 
Refs. [1,2, 3] for recent progress in this direction). 

For hadronic resonances, and in general systems involving 
two or more stable hadrons, the challenges are far greater and 
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further technological and formal developments are needed (see 
Refs. [4, 5, 6, 7] for recent reviews on the topic). In order 
to kinematically suppress multiparticle channels, many excited 
state calculations are performed using unphysically massive 
light quarks. Thus, it is desirable to devise a scheme for per¬ 
forming a controlled extrapolation to the physical mass. 

As a step towards developing such a program, we present the 
first extrapolation of a resonant scattering amplitude obtained 
from lattice QCD. Specifically, we analyze isovector, P -wave 
mi spectra in the elastic scattering region that have been deter¬ 
mined by the Hadron Spectrum Collaboration using dynamical 
quark masses corresponding to m„ = 236 MeV [8]. 

Lattice QCD uses a discrete and finite spacetime. Discretiza¬ 
tion provides a natural high energy regulator for QCD and if 
a fine enough spacing is used this introduces negligibly small 
effects in the spectrum. Working in a finite, periodic volume 
transforms the continuum of infinite volume scattering states 
into a discrete spectrum of states. The non-perturbative map¬ 
ping between finite and infinite volume observables was first 
derived in Refs. [9, 10] and is commonly referred to as the 
“Liischer method”. 

The mappings between finite and infinite volume amplitudes 
cannot be one-to-one due to two important facts. First, the re¬ 
duction of rotational symmetry from a continuous group to a 
discrete group (e.g., cubic) assures mixing between different 
partial waves. Second, having lost the notion of asymptotic 
states, finite volume states will necessarily be an admixture of 
different hadronic states with the same quantum numbers (e.g., 
nn and KK in the 7=1 channel). Many theoretical advances 
have guided the field. For example, several references have dis¬ 
cussed the feasibility of studying coupled-channel scattering in 
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a finite volume [11, 12, 13, 14] (see Refs. [15, 16] for the first 
application of this formalism to the study of nK, rjK) as well as 
three-body systems [17, 18, 19, 20], These methods become in¬ 
creasingly cumbersome when applied to highly energetic few- 
body systems, such as exotic or hybrid resonances [21, 22, 23], 
as well as the phenomenologically interesting charm and bot¬ 
tom decays (e.g., D —> nn/KK [13, 24]), where multiple few- 
body channels are open. 

In this work, we investigate one of the most studied low-lying 
resonances, the p [8, 25, 26, 27, 28, 29, 30, 31]. The p is an 
isotriplet with J pc = 1 , and it decays strongly to nn nearly 

100% of the time [32], Its mass, ~ 770 MeV, lies above the nn 
and 4n thresholds, and is less than half a width [T p ~ 145 MeV] 
away from the 6 n threshold. The coupling to these channels are 
experimentally observed to be negligible, which would suggest 
that the finite volume effects associated with these thresholds 
are suppressed. Further work is needed to confirm and quantify 
this suppression. 

To circumvent these subtleties, we perform an extrapolation 
to the physical point of the nn scattering phase shift computed 
at m n = 236 MeV [8]. At these quark masses, the 47r, 6n and 
KK thresholds lie well above the p resonance and can be safely 
ignored. To perform the extrapolation we use Unitarized Chi¬ 
ral Perturbation Theory (U\fPT) [33, 34, 35, 36, 37], which we 
summarize below. The parameters of U^PT at m„ = 236 MeV 
are chosen in order to reproduce the lattice QCD spectrum, 
and once this is done the pion mass is set to its experimen¬ 
tal value and a postdiction for the scattering phase shift is ob¬ 
tained. Although superficially the need to extrapolate may seem 
undesirable, the avoidance of thresholds makes this conjunction 
of a phenomenological effective field theory with the Liischer 
method a fruitful alternative to a determination of the phase 
shift at the physical point. 

U^PT was previously advocated in the literature as a tool to 
determine physical resonances from lattice QCD [38, 39, 40, 
41, 42, 43, 44, 45, 46], and it has been used in the study of 
the quark-mass dependence of the p mass [47] 1 . Instead of 
focusing on the pole of the resonant amplitude, which has been 
the main focus of previous chiral extrapolations, we fit the full 
resonant amplitude. Given the correlation between the energy- 
and quark-mass dependence of these amplitudes, we find that 
this is sufficient to obtain the quark-mass dependence of the 
amplitude and consequently its pole. 

In ref. [8], a total of 22 nn energy levels are obtained be¬ 
low the 4n/KK thresholds. Also determined are energy lev¬ 
els above these thresholds, and from them the KK phase shift 
and nn, KK inelasticity are obtained using the formalism first 
presented in [12, 13]. In this work, we analyze only the 
states in the elastic region. To relate these to an infinite vol¬ 
ume scattering amplitude, M(P), we use the generalization of 
Liischer’s formalism for two degenerate scalar particles in mov¬ 
ing frames [9, 10, 49, 50, 51] 

dzt[F- l (P,L) + M(P)] =0, (1) 


1 It also has been used to determine the low-energy coefficients (LECs) for 
heavy-light systems by studying the quark-mass dependence of the scattering 
phase shifts of weakly repulsive channels [48], 


where F(P,L) is a function that depends on the total four- 
momentum P and the spatial extent of the cubic volume L, and 
the determinant acts on the space of spherical harmonics (for 
an exact definition of these quantities see Ref. [50]). This ex¬ 
pression is exact up to exponentially suppressed corrections that 
scale as e~ m " L , which we can safely ignore given that m„L « 4.4 
for the lattice used [8] 2 . Because the two particles are degen¬ 
erate, odd and even partial waves do not couple, even when the 
system is in flight. Furthermore, in Ref. [8] it was shown that 
in the elastic region the { > 3 phase shifts are consistent with 
zero. Therefore, Eq. 1 effectively gives a one-to-one relation 
between the spectrum and the elastic (l, I) - (1,1) nn scatter¬ 
ing amplitude. For real values of the relative momentum, q, the 
inverse of the scattering amplitude is related to the scattering 
phase shift 6 in the standard way [50] 


q cot <5^ = 1 6nE* n Re 




( 2 ) 


where E* n = 2 ^jq 2 + m 2 is the total energy in the center of 
mass (c.m.) frame. 

We use SU(2) U^PT to obtain the nn amplitude. Just like 
standard ^PT [53, 54, 55, 56, 57, 58], U^PT allows one to 
evaluate observables analytically in a perturbative expansion 
defined by (m n /4nf n ) 2 , where f n = 92.2 MeV [32] is the de¬ 
cay constant of the n. At each order in the expansion, one can 
write the scattering amplitude as a function of a finite number of 
LECs. At leading-order (LO) in the expansion only two LECs 
appear (m 0 and /o). At next-to-leading order (NLO) four other 
LECs emerge ({ r i=l _ 4 ). See Appendix A for the Lagrangian as 
well as perturbative expressions for the pion mass, decay con¬ 
stant, and the pion-pion scattering amplitude. When perform¬ 
ing the fit to the lattice spectrum, we fix mo such that m n = 236 
MeV. Given that the decay constant has not been determined, 
/o is fixed to reproduce the experimental value of f n . 3 The 
l r . cannot be directly obtained from the physical values of the 
mass and decay constant, but can be accessed from the scatter¬ 
ing amplitude. For the l — 1 partial wave, only two linear com¬ 
binations of these are needed to describe the scattering phase 
shift (a i = -2t r x + t r 2 and an = fjj). As discussed below, we 
fix these parameters by performing a fit to the lattice spectrum. 
Although the Pi are quark-mass independent in principle, by 
ignoring higher-order corrections the LECs will absorb a mild 
quark-mass dependence. See Ref. [60] for a recent review and 
discussion in the context of standard^PT 4 . 

The distinguishing feature of UyPT is its use of a procedure 
commonly referred to in the literature as the Inverse Amplitude 
Method [33, 35, 36] to ensure that the scattering amplitude sat¬ 
isfies unitarity. Effectively, in U^PT ^-channel diagrams are 
summed in a geometric series using perturbation theory to all 


2 A subset of these exponential corrections has been determined for the nn 
states with £ = 0 [52] and ( = 1 [38] partial waves. 

Tor progress towards determining the decay constant of the ground state 
and excited states of the n using these lattices, we point the reader to Ref. [59]. 

4 In Ref. [47] it is argued that these effects might be large for UyPT and 
higher order corrections might be needed. In this work we ignored higher order 
corrections, and these will be incorporated in future studies. 
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Figure 1: Shown is the values of a t mQ previously determined for four dif¬ 
ferent values of quark masses (black circles) [65, 8]. The green band de¬ 
picts the fit to these masses using Eq. 4. The physical point is denoted by 
the dashed line. By fixing the resulting value of a t mQ to a t m^ ys ' we obtain 
4 2] = 0.1630(14) GeV' 1 . 


Figure 2: Shown is the I = \ irn phase shift obtained from the lattice QCD spec¬ 
trum determined at m n = 236 MeV as a function of the c.m. energy. The band 
corresponds to the SU(2) IJ^PT fit. The dashed line shows the 4 n threshold. 
We do not show two noisy energy levels. 


orders, while t- and w-channel diagrams are treated perturba- 
tively to a finite order in the expansion described above 5 . This 
procedure empirically extends the range of applicability of stan¬ 
dard *PT to c.m. energies on the order of 1.2 GeV. Further¬ 
more, unlike standard ^PT, UyPT has been shown to accu¬ 
rately describe low-lying resonances with a finite number of 
LECs [33, 35, 36], making it a desirable tool for the study of 
resonances from lattice QCD. By truncating the chiral expan¬ 
sion to NLO, one can write the unitarized scattering amplitude 
(see Appendix B for the derivation), 


WlufPT - Afl.O 


1 

W(lo - WInlo 


Af[,o, 


(3) 


where WIlo and WInlo are the LO and NLO %PT amplitudes 
detailed in Appendix A. 

To perform a chiral extrapolation we must determine the lat¬ 
tice spacing. We use two definitions of the lattice spacing. 
First, we use the O baryon mass, which has been determined 
to be = 0.2789(16) at these quark masses [8], By set¬ 

ting this equal to a,m p ^ ys ', where m p ^ yh ' = 1672.45(29) MeV 
is the mass of physical Q baryon, we obtain the lattice spac¬ 
ing = 0.1668(10) GeV -1 . Second, as shown in Fig. 1, we 
perform an extrapolation to the physical point of the lattice Q 
baryon mass using 

2 4 

mt ml 

m n (m n ) = m n .o + a—+ p—r (4) 

m n m n 

determined for four different values of e [0.14 - 0.33] [65, 

8], We find <ip ] = 0.1630(14) GeV -1 with a* 2 * /d.o.f. = 0.52. 
Assuming that aj 11 should coincide with a{ 2] , we perform all fits 
using both of these lattice spacings and any deviation of the re¬ 
sult is incorporated into the systematic error. All central values 
below are obtained using the mean value of a,. As shown be¬ 
low, this 2% error is the largest source of uncertainty in our final 


result. It is important to recognize that this systematic error is 
improvable. 

We determine the two unknown LECs by fitting the 22 en¬ 
ergy levels obtained at a single quark mass and spatial volume. 
In practice, we input the UyPT amplitude into Eq. 1 and com¬ 
pute the spectra for a given set of LECs, £ u -r pT ([a,l). By vary¬ 
ing these LECs we minimize the defined as 

* 2 ({<*/)) - X 8E M^\) C ■{ «£*({<*,-}) (5) 

j,k 

where 6Ej({ai}) = [.E*? 1 — £ , y^ PT ([a',])J, and {j, k) run over all 
22 energy levels. As with the energy levels themselves, the ele¬ 
ments of the covariance C matrix were provided by the Hadron 
Spectrum Collaboration [8]. The fit results in^ 2 /Vd. 0 .f. = 1-26 
for SU(2) U^PT and is shown in Fig. 2 compared to the lattice 
determined phase shifts. The LECs and correlations are found 
to be 

<*i(770 MeV) = 14.7(4)(2)(1) x 1(T 3 
<*2(770 MeV) = -28(6)(3) (°j) x 10 -3 

The first uncertainty is statistical, the second is the systematic 
due to the determination of the n mass and the anisotropy of the 
lattice 6 , and the third is an estimate of the systematic due to 
the determination of the lattice spacing. The symmetric matrix 
on the right of the coefficients denotes the statistical correlation 
between the two. By analytically continuing the scattering am¬ 
plitude to complex values of ,v = (E* n ) 2 we obtain a resonance 
pole on the unphysical sheet, corresponding to taking the neg¬ 
ative root when computing the c.m. momentum q*. At these 
quark masses, we find a p pole at E p = 782(2) - j 85(2) MeV 
with a width, T p = -2 Im(E p ) = 85(2) MeV. We observe good 
agreement with the result from the Hadron Spectrum Collabo¬ 
ration where the poles were determined using other parameter- 
izations of the scattering amplitude. This emphasizes the fact 


-0.98 

1 


( 6 ) 


5 We point the reader to Ref. [61, 47] for a rigorous derivation using dispel'- 

sive techniques [62, 63, 64]. The authors are not aware of such a derivation for 

inelastic processes, e.g., jiK —» t]K [15, 16]. 


6 The 7i mass was determined in lattice units to be a t m n = 0.03928(18). 
The anisotropy of that lattice is defined as £ = a s la t where a s and a r are the 
lattice spacings in the spatial and temporal extents. The anisotropy has been 
determined to be £ = 3.4534(61). 
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Figure 3: I = 1 7T7t phase shifts at three pion masses. In red we show the 
lattice-determined phase shifts, along with the SU(2) U^PT fit to the spectrum 
at m n = 236 MeV. The green band shows the extrapolation to the experimental 
pion mass. In blue we show the discrete points from the lattice calculation at 
m n - 391 MeV [26] and the extrapolation from the parameters determined from 
this 236 MeV fit. The extrapolated bands include both statistical and systematic 
errors discussed in the text. 


that the lattice QCD spectrum properly constrains the scattering 
phase shift independently of the parameterization chosen. 

The power of the U^PT amplitude is that it allows one to ex¬ 
trapolate these quantities as a function of pion mass. In Fig. 3 
we show the result of this exercise using the mean values of the 
coefficients in Eq. 6 and propagating both statistical and sys¬ 
tematic uncertainties. We show the postdiction for m n = 140 
MeV and m n = 391 MeV, where an earlier calculation also 
extracted the nn scattering amplitude containing the p reso¬ 
nance [26]. We emphasize that in Ref. [47] it is clearly ex¬ 
plained that UyPT is not expected to reliably describe lattice 
QCD results above m n ~ 300 - 350 MeV. Despite this formal 
constraint and the slight deviation at m n = 391 MeV from the 
lattice results, U^PT produces phase shifts that resemble both 
experimental and lattice determinations as a function of m n . 

In Fig. 4 we show a comparison of the results of the extrap¬ 
olation using SU(2) and SU(3) versions of IJyPT. Given that 
SU(3)-breaking effects are large, SU(3)^'P7’ has a poorer con¬ 
vergence than that the SU(2) counterpart. Therefore, we expect 
the SU(3) extrapolation to have a significantly larger systematic 
uncertainty. Assessing such systematic lies outside of the scope 
of the present work. 

In Fig. 5 we present our final result for the chiral extrapo¬ 
lation of the nn phase shift using SU(2) UyPT. The result in¬ 
cludes a propagation of statistical and systematic uncertainties. 
The largest uncertainty is due to the determination of the lat¬ 
tice spacing, where we aim to be conservative. Overall, we find 
good agreement with the experimental phase shift [66, 67] up to 
center of mass energies of 1.2 GeV, well above the An, 6n, KK 
and 8/r thresholds. By analytically continuing the amplitude 
into the complex plane, we find a postdiction of the p pole at 
the physical point E p = [755(2)(1)( 3 °) - | 129(3)(1)(’)] MeV. 

In order to compare with experimental determinations of the 
mass and width of the p, we must restrict out attention to those 
determinations which have used the model-independent defi¬ 
nitions nip = Re(£p) and I ' p = -21m (E p ). We contrast this 
with the standard procedure of quoting the mass and width 



Figure 4: Extrapolations of the phase shift determined at m n ~ 236 MeV 
(green) to the physical point done using SU(2) and SU(3) in blue and red re¬ 
spectively. The extrapolated bands include only statistical error. For an estimate 
of systematics and a comparison with experimental data, see Fig. 5. 


parameters appearing in the Breit-Wigner parametrization of 
the scattering amplitude (as is done in the Particle Data Group 
book [32]). Only in the very narrow width limit do these two 
definitions coincide. 

In Fig. 6 we show our determination of the p pole. For com¬ 
parison we show those obtained in Refs. [68, 62, 54, 69, 70, 71] 
by solving the Roy equation [72] and using experimental data as 
input. Since these results cover a large area, we highlight a dark 
point which encompasses all pole positions. Identifying this 
as an estimate of the overall systematic and statistical uncer¬ 
tainty, we find good agreement with our determination. We also 
show the pole position obtained in previous lattice QCD calcu¬ 
lations [65, 26, 8], including those where the p is stable. This 
plot serves as a nice illustration of the trajectory being taken by 
the p pole as a function of m n . For heavy quark masses, the p is 
stable and its pole lies on the real axis. As the quark mass de¬ 
creases, the p becomes unstable and acquires a non-zero width, 
sending the pole off the real axis. 

We compare the LECs determined here with those deter¬ 
mined in Refs. [46, 47, 37]: oa(770 MeV) x 10 3 6 [9,13] and 
a2(770 MeV) x 10 3 e [1,12]. We observe a qualitative discrep¬ 
ancy between our determination of a 2 and those determined in 
these references. This can be explained by two facts. First, as 
discussed in Ref. [34], the {(, I ) = (1,1) amplitude primarily 
depends on a\. Second, as mentioned above, the definition and 
value of these parameters depend on higher order corrections 
in the chiral expansion [47]. We suspect that by performing 
simultaneous fits of various channels while including higher 
order corrections one will see a convergence of these results. 
Implementing these techniques for channels including scalar 
resonances like the /o(500) would require using the modified 
Inverse Amplitude Method to have the correct analytic structure 
below threshold [61, 47], The implementation of this awaits the 
lattice QCD calculation of these channels using m n - 236 MeV. 

Final remarks: We present the first extrapolation of a reso¬ 
nant amplitude from lattice QCD. To perform the extrapolation 
we used UyPT, an effective field theory that at low-energies 
coincides with ^PT and at high-energies generates resonances 
dynamically. In this framework, resonances are manifested nat- 
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Figure 5: Shown is the extrapolation to the physical quark masses of the (£, /) = (1,1) nn scattering phase shift. This is plotted as a function of the c.m. energy 
( E * n ). The darker blue inner band includes only statistical uncertainty, while the lighter outer band also includes systematic uncertainties explained in the text. We 
see good agreement with the experimental phase shift shown as black circles [66] and green squares [67]. The dashed lines denote the 4 n, 6 n, KK and 8 n thresholds, 
which appear to play a negligible role. 
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demonstrate that by properly constraining the scattering ampli¬ 
tude at a value of the pion mass where fewer channels are kine¬ 
matically open, one can perform an extrapolation to the physi¬ 
cal point. 

These methods may be applied to obtain a wide range of 
hadron scattering amplitudes that are presently being extracted 
from lattice QCD, in both the light and heavy quark sec¬ 
tors [26, 73, 74, 15, 16, 28, 75, 76]. It is hoped that these 
concepts could be extended and applied to scattering processes 
containing highly excited and exotic resonances to gain deeper 
understanding of QCD and the excited spectrum of hadrons. 


Figure 6: We compare our determination of the p pole [red square] with previ¬ 
ous lattice calculations [blue circles] performed using unphysically heavy quark 
masses where the p is stable [65] as well as unstable [26, 8]. This is compared 
with pole determinations obtained from solutions to the Roy equation [72] con¬ 
strained from experimental data [gray diamonds] [68, 62, 54, 69,70, 71]. Using 
these we highlight a black diamond whose uncertainty is defined to include all 
determinations from these references up to one standard deviation. 


urally as singularities in amplitudes. We observe that this ef¬ 
fective held theory does a remarkable job in describing the re¬ 
cent results of the Hadron Spectrum Collaboration. Using the 
lattice QCD spectrum to constrain the LECs of the theory, we 
find good agreement with the experimentally measured nn scat¬ 
tering phase shift up to energies above the 87r and KK thresh¬ 
olds illustrating the significance of this result. We observe the 
extrapolated amplitude to have a pole, corresponding to the p 
meson, which agrees with previous determinations using dis¬ 
persive analysis of experimental nn scattering data. 

It is desirable to study more complex systems such as highly 
energetic exotic hadrons (e.g., the 7Ti(1400) resonance) or heavy 
meson weak decays (e.g., D —> nn/KK [13, 24]), however it is 
not yet clear when a finite volume formalism rigorously accom¬ 
modating all open multiparticle channels will be available. We 


Appendix A. Chiral Lagrangian and Scattering Ampli¬ 
tude 

Here we present the key results of SU(2) ^P'f as derived in 
Ref. [58]. The relevant terms of the leading order (LO) and 
next-to-leading order (NLO) terms of the chiral Lagrangian (in 
the isospin limit m u = mp, 

f 2 nr f 2 

-Clo = J ~j Tr^f/^t/ 1 ) + Tr(f/ f + U ) 

-Cnlo = j [Tr^f/^l/)] 2 + | [Tr (d M Ud v U f )] [Tr (d fl Ud v U t 

+ [Tr (U 1 + U)] 2 + Tr (d 2 l 7 f - 3 2 (/) + ..., 

(A.l) 

are written in terms of the parameters fo (related to the pion de¬ 
cay constant) and mo (related to the pion mass), and the matrix 
of pion fields, 


U = exp 



y[2n 



(A.2) 
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Divergences associated with loops with LO vertices are re¬ 
moved by renormalizing the 7; LECs from the NLO Lagrangian 
and physical quantities depend on the renormalized LECs 
We use jj = 770 MeV in this work. At this order in the chiral 
expansion, it is convenient to introduce //-independent expres¬ 
sions for the LECs, that depend on the value of mo. 




7 i 

32 7T 2 


\(i + In , 


(A.3) 


where y x = ~, y 2 = y 3 = and y 4 = 2. 

We use the standard NLO expressions [56, 57, 58] for the 
physical pion mass and decay constant. 


1 ml 


32* 2 / 0 2 


{3 (mo) + ■■■ 


fn = fo 


1 + 


16tt 2 f- 




(A.4) 


(A.5) 


to solve for m () and fo perturbatively. To fix mo we use the value 
of m n that has been determined on the lattice, m ] f a '. Since f„ 
has not been determined for these lattices, we resort to fixing f„ 
using the experimental value, . This approximation forces 
us to use two different values of m„ in our fits. More explicitly, 
for mo we use. 


ml * (m l r f 


1 + 

1 + 


1 (mjf u ) 2 

32 ^-—fT 

1 (m'f 1 ) 2 

32^ 2 (/“ p ) 2 


+ ■ ■ ■ 

h (mJT) + ■ ■ ■ 


(A.6) 
(A.7) 


were the ellipses denote corrections that appear at higher orders 
in the chiral expansion. Similarly, for fo. 


fo * /f 



1 (m“ p ) 2 

i6^ 2 ( /r p ) 2 


4(w“ p ) +... 


(A.8) 


The amplitudes depend on 1 // 2 , which we write here perturba¬ 
tively 


where 



(A. 11) 

Note that in Eq. A. 10 we have implemented the perturbative 
7 1 

expressions for ;n[j and — described above. In Eq. A. 10 and 
fo 

Eq. A. 11 we use the notation m n = m^ lt and f„ = /™ p '. The 
amplitude ( 5 , t, u) can then be projected into a partial wave (' 
using. 


1 r 1 

Mf = - J dzPeiz) Jl(s,t(s,z),u(s,z)) (A. 12) 

where z = cost? and 0 is the ,v-cflannel c.m. frame scattering 
angle. In this work, we also project onto the 7=1 channel, 

M l (s, t , u ) = Al(f, s, u ) - A4(u, t , s). (A.13) 

One can show that the only linear combinations of LECs con¬ 
tributing to the isotriplet scattering amplitude are a\ = -27] +{ r 2 
and which are the ones determined in this work. 


Appendix B. The Inverse Amplitude Method 

Although UyPT has been extensively discussed in the lit¬ 
erature, here we sketch the derivation of Eq. 3 presented in 
Ref. [36] in an effort to make this article more self-contained. 
The basic idea, as already mentioned above, is to assure that 
unitarity is satisfied exactly at each order in the chiral expan¬ 
sion. We begin by giving the standard relation between the 
S -matrix and the partial-wave projected scattering amplitude, 
M, 

5 = 1+ 2/ixM, (B.l) 

where cr = qj l 6 nE* n . Unitarity enforces 


f 2 o U ? V ) 2 


2 (m exp ) 2 

^ v m 7 T ) 75 / exp. 


16* 2 ( /; xp ) 2 


UmT) + • ■ ■ 


(A.9) 


The scattering amplitude prior to partial-wave projection, 
2 H(s, t, 11 ), can be written as 


Jl L o(s, t, u) = 


^1nlo(A t, u) = 


s - mz 


/2 

s - ml (m“ p ) 2 - 


Um7 P ) - 


fn 2 87T 2 / 2 * 327r 2 /4 

{3(s 2 - ml)J(s) + [f(f - u) - 2 mil + 4w 2 m - 2 ml | /(f) 

+ | u(u — t) — 2 m 2 n u + 4m 2 n t - 2n?^J J(u )J 

{2(7 1 (»0-4/3)( i -2m 2 ) 2 

+( 62 (m n ) - 5/6) [s 2 + (f - m) 2 ] - 12 mis + 15«4}, 

(A. 10) 


96^ 2 /4 


Im(Al) = o-|M| 2 , 


(B.2) 


which is the familiar Optical Theorem. This condition can be 
rewritten as 

ImfAT 1 ) = — cr, (B.3) 


which leads us to 


M = (Ref AT 1 ) - kr )~ 1 . 


(B.4) 


If A1 is evaluated perturbatively as detailed in Appendix A, 
At = Mu) + AInlo + ..., we can expand its inverse to find. 


AT 1 = AT 1 


1 


L0 1 + A1 L q AInlo + ■ • • 


- ^lo (l “ A^lo AInlo + •■■)■ 


(B.5) 


Since AIlo is real, 

Re(AT') = Mlo (l - A4 L o Re(At NL o) + •■■), (B.6) 
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which we insert into Eq. B.4 to find. 


^LO (l ~ ^LO ^ e (^NLo) +•■•)- icr 

~ AJlo-;-^—WIlO- (B.7) 

/VIlo - RefA^NLo) _ icrM ^ 0 

Finally, let us return to Eq. B.2 and enumerate the unitarity 
constraints order by order, 

LO : Im(M L0 ) = 0 

NLO : Im(M NL0 ) = cr|M LO | 2 = <rM 2 0 . (B.8) 

Thus, putting Eq. B.8 into Eq. B.7, we reproduce Eq. 3. 
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